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Abstract 

Background and Aims: Ticks are highly susceptible to global environmental and socio-economical changes. Several tick- 
borne pathogens have been reported in new geographical regions while new species, strains or genetic variants of tick- 
borne microorganisms are continually being detected. However, tick-borne pathogens are still poorly understood, and it is 
estimated that half of all human tick-borne disease has an unknown origin. Therefore in order to prevent these diseases, 
more effort is required to identify unknown or unexpected tick-borne pathogens. Ixodes ricinus is the vector for a broad 
range of bacterial pathogens and the most prevalent tick in Europe. The aim of the present study was to evaluate the 
capability of Next Generation Sequencing (NGS) to extend the inventory of pathogenic bacteria carried by this species of 
tick in France. 

Methods: RNA and DNA were extracted from 1450 /. ricinus questing nymphs collected by flagging in Alsace, France. RNA 
was pooled and used for NGS. Following de novo assembly, bacterial contigs were assigned to the closest known taxonomy. 
DNA was used for real time PCR to confirm taxonomic species assignment of NGS-derived contigs for the doubtful cases, 
and for determination of prevalence. 

Results:\Ne have generated a global in-depth picture of tick-borne bacteria. We identified RNA from the main pathogenic 
bacterial species known to be transmitted by /. ricinus. In addition we also identified unanticipated bacterial species for 
which we have estimated the prevalence within those ticks inhabiting the studied areas. 

Conclusions: The data obtained from this study has proven that NGS has an enormous potential to detect the unexpected 
and provides the means to monitor pathogen occurrence. 
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Introduction 

Ticks are highly susceptible to global environmental and socio- 
economical changes [1], such that the incidence of tick-borne 
disease is increasing worldwide [2], and that some tick-borne 
pathogens have been reported in new geographical regions [3,4]. 
Ixodes ricinus is the most common European tick species, and is the 
vector for a broad range of bacterial pathogens [2,5]. 

Due to recent advances in molecular biology, new species, 
strains or genetic variants of microorganisms are continually being 
detected in ticks, resulting in an ever-increasing list of pathogens 
capable of infecting livestock and humans. In Europeans, the most 
prevalent tick-borne human disease is Lyme borreliosis, caused by 
any of four different Borrelia species; B. burgdorferi sensu stricto, B. 
afzelii, B. garini, and B. spielmanii [6], with over 50 000 new cases 
each year in Europe [7]. Lyme disease has a high rate of diagnosis, 
mainly due to its distinctive symptoms and the existence of 
effective serological tests. Although, in addition to these typical 



cases, patients with tick bites can also be infected by many other 
bacteria including Anaplasma phagocytophilum, Rickettsiae spp., 
"Candidatus Neoehrlichia mikurensis", Ehrlichia spp., Bartonella 
spp. and Francisella spp. Some of these species have been identified 
in ticks decades prior to their association with human disease, 
whereas others have only been discovered very recently. Therefore 
we can confidendy assume that many tick-associated pathogens 
are still unknown [3,4,8]. In cases of co-infection, symptoms 
induced by these microorganisms are often non-specific and can 
be aggravated by pre-existing or concurrent diseases [9] making 
diagnosis exceedingly difficult. It is estimated that up to half of 
human tick-borne disease has an unknown microbiological origin 
[10]. Consequendy, in order to identify novel or unexpected tick- 
borne pathogens in their natural hosts, concerted efforts are 
required to develop efficient detection methods. 

Up until now, most studies aiming to identify tick-borne 
pathogens in their vectors have only been able to assess minimal 
numbers of pathogens at a time, which is partly due to 
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technological limitations. Indeed, complete screens of microor- 
ganisms in natural host populations (both reservoirs and vectors) 
were out of reach using standard laboratory procedures. Identi- 
fication of microorganisms in biological samples has been 
dominated by the use of culture-dependent methods and by 
conventional molecular approaches. These methodologies suffer 
from major limitations: (i) isolation in cell culture or in synthetic 
medium is difficult for many microorganisms and (ii) molecular 
identification was mostly based on the use of specific primers 
combined with real-time PCR, which are only capable of detecting 
a selected number of species. However, the recent rapid 
development of Next Generation Sequencing (NGS) methods 
combined with bioinformatics has revolutionized the research field 
of epidemiology and diagnosis of infectious diseases, and has 
effectively overcome those limitations. NGS could be effective 
tools in diagnosing pathogens for medical and veterinary purposes 
[11-18] and these techniques has recently been successfully used 
to identify the bacterial communities associated with I. ricinus ticks 
[19,20]. In these studies, the authors used the amplification and 
the sequencing of the hyper-variable region V6 of the 16S rRNA 
encoding genes. If this gene is well adapted to describe bacterial 
communities, the analysis of its sequence does not allow identifying 
bacteria at the species level. Ticks carry many bacterial species, 
including endogenous ones (symbionts and commensals). In ticks, 
symbionts are often phylogenetically close to pathogens and 16S 
rRNA diversity does not allow, for instance, distinguishing 
symbiotic Rickettsia species from pathogenic ones. Besides this 
problem, within a bacteria genus, many species exist. This is the 
case for Bartonella sp., for which some species are known to be 
pathogenic for humans (or animals) and other have never been 
associated to diseases. Thus identifying them at the species level is 
very important. This is also the case for the bacteria belonging to 
the complex of Borrelia burgdorferi sensu lato, responsible of Lyme 
diseases, for which the analysis ofl6SrRNA diversity is not efficient 
to identify them at the species level. 

In this context, we undertook a study to produce an inventory of 
predicted and/ or unexpected pathogenic bacteria carried by ticks, 
in a geographical region with abundant ticks and a concomitant 
high risk of disease transmission. To overcome the limitation of 
deep sequencing of 16S RNA encoding genes, we sequenced the 
whole transcriptome of ticks and their resident microbiota. 

Materials and Methods 

Tick collection and extract preparation 

No specific permits were required for the described field studies. 
These field studies did not involve endangered or protected 
species. A total of 1450 questing nymphs and 62 females of the I. 
ricinus tick were collected by flagging from three wooded habitats 
in Eastern France (Alsace Department, [Murbach (47°55'05"N, 
7°8'46"E), Hohbuhl (48°27'33"N, 7°17'22"E) and Wasselonne 
(48°38'09"N, 7°21'45"E)]). All collected ticks were washed as 
previously described [21], nymphs were pooled into groups of 25 
individuals (58 pools in total) and adults into groups of 2 
individuals (31 pools in total) and crushed in 300 ul of Dulbecco's 
MEM (DMEM) medium supplemented by 10% fetal bovine 
serum. A pool of 15 /. ricinus nymphs from our pathogen-free 
colony was equivalently treated and used as a reference. This 
control colony originated from female ticks collected in Murbach 
(Alsace, France) and was bred as previously described [22]. The 
pathogen-free colony was obtained as follows; females from the 
field were engorged on rabbits and allowed to lay eggs. Female tick 
DNA samples were extracted and different PCRs were performed 
to test for the presence of Borrelia spp., Bartonella spp., Anaplasma 



spp., Rickettsia spp., Francisella spp. and Coxiella spp. Only larvae 
from "pathogen-free" female ticks were conserved and maintained 
in our colony. 

DNA and RNA extraction 

Crushed tick pools were divided into two equal samples. Total 
RNA, used to identify viable and replicating microorganisms, and 
total DNA, for both pathogen taxonomic confirmation and 
prevalence studies, were then extracted from the nymphs using 
the Nucleospin RNA II kit (Macherey-Nagel, Duren, Germany), 
and the Wizard genomic DNA purification kit (Promega, 
Madison, WI, USA) following manufacturers' instructions. 

High throughput sequencing and data analysis 

High throughput sequencing was only performed on nymphs, 
and all RNA samples were pooled before sequencing. RNA was 
first retro-transcribed to cDNA, then randomly amplified using the 
Multiple Displacement Amplification (MDA) protocol with phi29 
polymerase and random hexamers as described [17]. Library 
preparations and sequencing with an Illumina HiSeq2000 
sequencer were outsourced to DNA Vision (Belgium). Wild and 
pathogen-free samples were sequenced to a depth of 1 00 Mio and 
62 Mio paired-end reads of 101 bp respectively. Raw sequence 
reads were trimmed according to their quality score. At the time of 
analysis, there was no published reference genome for /. ricinus 
available in public databases, so sequences from the ticks 
themselves, or from symbiotic or commensal bacteria naturally 
found in ticks, were removed from the analysis by subtracting 
sequences derived from the pathogen-free reference sample using 
the SOAP2 aligner tool. Finally, de novo assembly was performed 
on all remaining reads (7.7 Mio), producing 174 841 contigs. 
Taxonomic assignation of these contigs was achieved via 
successive sequence alignment using nucleotide and protein 
databases and the BLAST algorithm. Contigs were assigned the 
closest known taxonomy according to their identity percentage, 
and distant alignments were not considered. Of the assigned reads, 
12.22% of the cDNA derived sequences were of bacterial origin. 

Specific Real-Time PCR for taxonomic species assignment 
of NGS-derived contigs confirmation 

Real-time PCR on DNA samples was performed to confirm 
taxonomic species assignment of NGS-derived contigs for those 
cases in doubt: Francisella tularensis, Coxiella burnettii, Ehrlichia canis, 
Rickettsia felts and R. canadensis. All DNA samples were pooled as for 
RNA analysis. Primers and probes used for amplification are 
described in Table 1. 

DNA pre-amplification was performed using the TaqMan 
Pre Amp Master Mix kit (Applied Biosystems). Primers were 
pooled at a final concentration of 0.2X with equal primer 
concentrations (20 uM each). The reaction was performed in a 
final volume of 5 (0.1 containing 2.5 ul TaqMan PreAmp Master 
Mix (2X), 1.2 ul of the pooled primer mix (0.2X) and 1.3 ul of 
DNA, with one cycle at 95°C for 10 min, 14 cycles at 95°C for 
15 sec and 4 min at 60°C. After the pre-amplification PCR, 
products were diluted 1:20 and stored at — 20°C prior to use in 
qPCR. All fluorogenic probes were synthesized with a 6-carboxy- 
fluorescein (FAM) reporter molecule attached to the 5' end and a 
Black Hole Quencher 1 (BHQ1) attached to the 3' end. Real-time 
Taqman PCR assays were performed in a final volume of 12 ul 
using the LightCycler® 480 Probe Master mix (Roche Applied 
Science, Germany) at a IX final concentration, with primers and 
probes at 200 nM and pre-amplified DNA. Thermal cycling 
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Table 1. Primers and probes used for pathogen confirmation or prevalence estimation in ticks. 



Pathogen 



Target Primers and Probe Sequence (5' -> 3') 



Amplicon size (bp) Reference 



Borrelia miyamotoi 



Coxiella burnetti 



Ehrlichia canis 



Francisella tuiarensis 



Rickettsia feiis 



Rickettsia canadensis 



glpQ 



idc 



dsb 



tu!4 



orfB 



ompA 



Candidatus Neoehrlichia mikurensis groEL 



CACGACCCAGAAATTGACACA 

GTGTGAAGTCAGTGGCGTAAT 

TCGTCCGTTTTCTCTAGCTCGATTGGG 

AGGCCCGTCCGTTATTTTACG 

CGGAAAATCACCATATTCACCTT 

TTCAGGCGTTTTGACCGGGCTTGGC 

AATACTTGGTGAGTCTTCACTCA 

GTTGCTTGTAATGTAGTGCTGC 

AAGTTGCCCAAGCAGCACTAGCTGTAC 

ACCCACAAGGAAGTGTAAGATTA 

GTAATTGGGAAGCTTGTATCATG 

AATGGCAGGCTCCAGAAGGTTCTAAGT 

CCCTTTTCGTAACGCTTTGCT 

GGGCTAAACCAGGGAAACCT 

TGTTCCGGTTTTAACGGCAGATACCCA 

AGTGGGACATTGGGTGTTGC 

GGCAACACCTTGAGCATTATC 

CTCCTGCCTGCGTTATACCATGCCA 

AGAGACATCATTCGCATTTTGGA 

TTCCGGTGTACCATAAGGCTT 

AGATGCTGTTGGATGTACTGCTGGACC 



94 



74 



110 



96 



This study 



Adapted from [47] 



Adapted from [30] 



Adapted from [48] 



[49] 



This study 



This study 



doi:1 0.1 371 /journal.pone.0081 439.t001 



conditions were as follows: 95°C for 5 min, 45 cycles at 95°C for 
10 s and 60°C for 15 s and then 40°C for 10 s. 



from the positive pools were positive and prevalence if 
individuals from the positive pools were positive. 



1/2 



Specific Real-Time PCR for prevalence of unexpected 
pathogens in ticks 

Real-time PCR was performed to determine the prevalence of 
three unexpected tick-borne bacteria DNA in ticks: "Candidatus 
Neoehrlichia mikurensis", Borrelia miyamotoi, and R. fills. We 
assessed 58 pools of 25 nymphs according to the collection location 
(corresponding to a total of 1450 nymphs) and 31 pools of 2 
females (corresponding to a total of 62 females). The primers for 
amplification of "Candidatus Neoehrlichia mikurensis" and Borrelia 
miyamotoi were specifically designed for this study and are 
presented in Table 1 . Real-time PCR was conducted as described 
above. 

Statistical analysis 

For the pooled nymphs, we estimated prevalence at the 
individual level using the exact method of Hauck, assuming 
perfect sensitivity and specificity of our pathogen detection 
methods [23]. Hauck noted a one-to-one relationship between 
individual prevalence levels, 7t, and the prevalence of positive 
pools, P. A point estimate for the prevalence rate can thus be 
obtained from the pool positive rate by K = 1-(1-P) 1/A where kk the 
number of nymphs or adults per pool. Exact 95% confidence 
intervals were then obtained by assuming a binomial distribution 
for the number of positive pools [24] . For the adults, as each pool 
only contained two individuals, we calculated the individual 
prevalence intervals corresponding to prevalence if 2/2 individuals 



Results 

Inventory of medically important bacterial RNA content 
in ticks 

Bacterial content was organized into genera or putative species 
by comparison with the closest homolog (Figure 1). Species were 
only allocated unambiguous assignment to one specific species 
when either the NGS-derived contig identity percentage was 
&95% compared to a specific species sequence (and lower when 
compared to other species); or when NGS-derived contigs 
matched sequences from one single bacterial species. Following 
these criteria, all bacterial species unambiguously identified are 
listed in Table 2. 

Global analysis of genera mainly associated with medically 
important tick-borne bacteria: Anaplasma spp., Borrelia spp., Coxiella 
spp., Ehrlichia spp., Francisella spp. and Rickettsia spp., indicated that 
among the alpha-proteobacteria family, Rickettsia spp. were the 
most highly represented. Species identified included R. helvetica, a 
member of the Spotted Fever Group of Rickettsiae, predicted to be 
associated with human illness in many European countries [25]; R. 
canadensis isolated from ticks belonging to Haemaphy sails spp. in 
Canada [26], never previously identified in Europe, and involved, 
on serological evidence of infection, in Rocky Mountain Spotted 
Fever-like disease and cerebral vasculitis in US patients; and R. 
fills, responsible for Spotted Fever in humans and whose only 
known vector is the flea [27]. As the number of reads was low for 
R. canadensis (16) and R. fells (180), and because the presence of 
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Table 2. List of bacteria unambiguously identified at the species level in ticks, including the number of NGS reads for each species, 
the corresponding list of identified genes allowing taxonomic affiliation (with% identity) to one specific bacteria, and the closest 
bacterial species to the unambiguously identified species (with% identity). 



Unambiguously Identified bacterial species (read number) 


Gene hit (Accession number) -% id 


Closest species -% id 


Anaplasma phagocytophilum (5234) 


16S RNA (JN990105) - 100% 


A. ovis - 99% 




23S RNA (AF416766) - 100% 


A. marginale - 96% 


Bartonella henselae (33) 


dnaB (BH0270) - 100% 


fi. tribocorum - 94% 




dnaJ (BH00660) - 100% 


R quintana - 91% 




trkA (BH05860) - 100% 


B. quintana - 93% 




Phage related protein (BH09270) - 100% 


B. grahamii - 84% 




Hypothetical protein (BH022330) - 99% 


R tribocorum - 85% 




Transketolase - (BH15090) - 99% 


B. quintana - 87% 


Bartonella grahamii (39) 


Phage tail protein encoding gene (ACS5068.1) - 90% 


None 




Adhesin (ACS50528.1) - 75% 


None 


Borrelia afzelii (1057) 


Hypothetical protein (CP000395) 


B. garinii - 96% 


Borrelia garinii (723) 


Hypothetical protein (CP001205) - 100% 


fi. afzelii - 98% 




Hypothetical protein(CP003151) -100% 


fi. afzelii - 98% 


Borrelia burgdorferi s.s. (475) 


23S RNA (BORG23S) - 100% 


fi. garinii - 99% 


Borrelia miyamotoi (2499) 


16S RNA (JF951385) - 100% 


Other Borrelia sp. <98% 


Candidatus "Neoehrlichia mikurensis" (2010) 


16S RNA (JQ359061) - 100% 


A. phagocytophilum <97% 


Ehrlichia canis* (660) 


PhoT family ABC transporter (CP000107) - 75% 


None 




23SrRNA - 96% 


£ chaffensis - 96% 


Rickettsia canadensis (16) 


MraW (CP000409) - 96% 


R. fells - 93% 




TraD (CP000409) - 97% 


ft amblyommii - 85% 




16S RNA methyltransferase: (CP003304) - 97% 


ft montanensis - 88% 


Rickettsia fells (180) 


Plasmid pRF (GQ329881) - 98% 


ft australis - 91% 




Unknown protein (CP000053) - 95% 


ft australis - 90% 




AmpC (CP000053.1 ) - 99% 


ft amblyoma - 97% 




Sea 13 (CP000053) - 96% 


ft slovaca - 91% 


Rickettsia helvetica (2411) 


Seal (AY355363) - 100% 


ft fells - 91% 




ompB (JN036412.1) - 100% 


ft asiatica - 94% 




groEL (DQ442911) - 99% 


ft canadensis - 97% 




sca4 (FJ358501.1) - 100% 


ft tamurae - 94% 



doi:1 0.1 371 /journal.pone.0081 439.t002 



both pathogens in I. ricinus was not expected, we performed real- 
time PCR on DNA extracts corresponding to a single pool of 1450 
nymphs, which confirmed the detection of R.felis DNA, but not R. 
canadensis indicating that for this bacteria the RNA-seq results is 
below the threshold of detection. We can thus consider that this 
bacteria is not present in our samples 

Among the Anaplasmataceae family, we identified sequences 
belonging to A. phagocytophilum, the agent of human granulo- 
cytosis [25] and "candidatus Neoehrlichia mikurensis" a new tick- 
borne pathogen discovered in 2004. The first human case of 
"candidatus N. mikurensis" infection was reported in a Swedish 
patient in 2010, but has never been linked to human disease in 
France thus far [28]. The presence of DNA from this bacterium 
was confirmed by real-time PCR in the pool of 1450 nymphs. In 
2012, we had previously identified this bacterium for the first time 
in France, in bank voles from North Eastern France [29] . We also 
identified species belonging to the Ehrlichia genus. NGS sequence 
analysis combined with specific real-time PCR confirmation 
facilitated the detection of E. canis, the agent of canine 



ehrlichiosis, known to be transmitted by the dark tick Rhipiceph- 
alus sanguineus [30] (Table 2). 

Among the Bartonella spp. sequences, the cat-specific B. henselae 
was clearly identified, as well as sequences related to the rodent- 
specific B. grahamii (Table 2). Both of these species are known to be 
mainly flea-borne zoonotic human pathogens but have also been 
previously detected in European ticks [31-33]. 

We also detected RNA sequences related to gamma-proteo- 
bacteria belonging to Coxiella spp. and Francisella spp. For these 
genera, the low number of reads prevented the identification of 
these bacteria at the species level. We failed to detect Coxiella 
burnetii or Francisella tularensis DNA within ticks by real-time PCR, 
these two species being the most common tick-borne pathogens of 
each genera. 

As expected for ticks collected in Alsace, a French region with 
high rates of Lyme disease, we also identified B. afzelii, B. garinii, 
and B. burgdorferi s.s., the three human pathogenic species known to 
circulate in France [34]. Interestingly, we also detected sequences 
from B. miyamotoi, which is responsible for relapsing fever in 
Russia, and which was recently detected in the USA [35]. The 
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presence of B. miyamotoi DNA was confirmed by real-time PCR. 
Even though sequences from this bacteria have previously been 
amplified from I. ricinus in Sweden, France and Slovakia [21,36], it 
has never thus far been associated with disease in Europe. 

Prevalence of unexpected human pathogen DNA in ticks 

Considering their potential importance to human health and 
disease, we estimated the prevalence of "candidatus Neoehrlichia 
mikurensis", B. miyamotoi, and R.felis microorganisms. Among the 
58 tested pools of nymph DNA (each pool containing 25 nymphs), 
19 were positive for "candidatus Neoehrlichia mikurensis" and 25 
for B. miyamotoi. By estimating the individual prevalence of nymphs 
using the exact method of Hauck, we found a prevalence of 2.2% 
(1.4—3.3) for Borrelia miyamotoi and 1.6% (0.9-2.5) for "candidatus 
Neoehrlichia mikurensis". For adults, we analyzed 31 pools, each 
containing two adults. We identified two positive pools for 
"candidatus Neoehrlichia mikurensis" which represented an indi- 
vidual prevalence of between 3.2% and 6.4%. All tested adult 
pools were negative for B. miyamotoi. For R.felis, we identified one 
positive pool of nymphs among the 58 tested (estimated prevalence 
of 0.09%), but none within the adult pools. 

Discussion 

We have used a whole transcriptome analysis of ticks RNA to 
generate a global picture of tick-borne bacteria and validate an 
NGS strategy, by identifying bacteria at the RNA level. It has been 
recently shown that metagenomic 16S rDNA Illumina tags are a 
powerful alternative to 16S DNA amplicon sequencing [3 7]. This 
paper extends this approach by using wide transcriptome analysis, 
including but not limited to 16S and 23 S RNAs. We have 
identified bacterial species known to be transmitted by I. ricinus 
and causing disease in France or Europe: Rickettsia helvetica; 
Anaplasma phagocytophilum; Borrelia garinii, B. afzelii, B. burgdorferi s.s. 
and "candidatus Neoehrlichia mikurensis". Considering "candidatus 
Neoehrlichia mikurensis", has never been identified in ticks from 
France, we calculated an estimated prevalence of 1.6% in nymphs 
and from between 3.2 to 6.4% in adults, rates similar to those 
found elsewhere in Europe [3]. "candidatus Neoehrlichia mikur- 
ensis" infection has also been detected in humans from north- 
eastern China, and field surveys showed that 1.6% of ticks and 
3.8% of rodents collected from patient residences were also 
infected [38,39]. Although corresponding human or animal 
disease has not yet been reported in France, the prevalence of 
positive rodents and ticks is similar to levels found in China. In 
addition, the infected French bank vole (Myodes glareolus) and ticks 
are located in close proximity to human dwellings, suggesting that 
"candidatus Neoehrlichia mikurensis" is very likely to pose risks to 
public health in France. However, to date, the diagnosis of 
"candidatus Neoehrlichia mikurensis" infection has solely relied on 
PCR amplification of bacterial DNA. Therefore, the absence of 
serological tests combined with physicians' lack of information 
relating to these bacteria, makes diagnosis particularly difficult. 
Thorough surveillance and improved diagnostic tools are required 
to provide more insight into the relevance of Candidatus N. 
mikurensis to public health. 

In addition to the expected species, our study generated 
evidence of other unanticipated bacterial species. For example, 
B. miyamotoii, transmitted by Ixodes spp., and causing recurrent 
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fever, has not previously been associated with disease in Europe 
[35]. B. miyamotoi was isolated for the first time in Japan in 1995 
from I. persulcatus ticks as well as from the blood of Apodemus argentus 
mice [40,41]. In Europe, B. miyamotoi was detected in ticks from 
Sweden, Germany, and very recently, from Estonia [36,42,43]. In 
addition to A. argentus mice, it has been shown that other rodents 
such as white-footed mice (Peromyscus leucopus) may serve as host 
reservoirs for B. miyamotoi [44]. In our study B. miyamotoi was 
estimated to occur at a mean prevalence rate of 2.2% in nymphs, a 
rate close to those obtained from different European countries, 
which fluctuates from 0.4% in Estonia to 3.5% in Germany 
[35,36,43,45], Some symptoms caused by this particular Borrelia 
sp., such as recurring high fever, are very different from the skin 
rash and joint pain typical of Lyme disease. Identifying the RNA 
of this bacterium in French I. ricinus for the second time strongly 
suggests that the disease could exist in this country. However, in 
the absence of specific tests, it is impossible to know whether the 
disease is indeed present in France (or in Europe), and if it has the 
potential to cause serious long-term damage, as can result from 
untreated Lyme disease. 

We also detected known pathogens which are either normally 
transmitted by other tick species, or by vectors other than ticks, 
such as for the B. henselae, B. grahamii and R.felis microorganisms. B. 
henselae and R. felis are typically transmitted from cat to cat via 
fleas, with human contamination arising from cat or flea bites. 
Establishing ticks as vectors of B. henselae has been fiercely debated 
for many years, even given the recent demonstration of the 
competence of I. ricinus to host B. henselae [46]. To our knowledge, 
our data present the first characterization of R. felis RNA in ticks. 
However even though we have confirmed the presence of RNA 
and DNA from R.felis in these ticks by real time-PCR with specific 
primers, the detection of R.felis DNA in only one pool of nymphs 
strongly suggests, if it exists, an anecdotic role for /. ricinus in R. felis 
circulation. 

Detecting RNA of unexpected pathogens in Ixodes ticks indicates 
that these microorganisms are viable and actively replicating, 
however, it does not prove that transmission by I. ricinus actually 
occurs. Competence studies of /. ricinus for all detected pathogens 
should be experimentally clarified in the future before reclassifying 
them as tick-borne pathogens. Although further research is thus 
required to confirm the role of I. ricinus in transmitting these 
pathogens, the data obtained with NGS has proven that this 
method has enormous potential to detect the unexpected, and 
provides an excellent means for surveying pathogen occurrence in 
disease vectors from particular localities. 

Acknowledgments 

We thank the group Tiques et Maladies a Tiques group (TMT) of Reseau 
Ecologie des Interactions Durables for stimulating discussion and support, 
the DIM He de France, the CoVetLab partner institutes, the ANSES, the 
Animal Health department of the INRA, the INRA Metaprogramme 
MEM (Patho-ID) and EU FP7-261504 EDENext for their support. 

Author Contributions 

Conceived and designed the experiments: MVT ME. Performed the 
experiments: LM ED JC CH. Analyzed the data: MVT SM SB JC CH 
ME. Wrote the paper: MVTJC ME. 



2. Parola P, Raoult D (2001) Ticks and tiekborne bacterial diseases in humans: an 
emerging infectious threat. Clin Infect Dis 32: 897-928. 



PLOS ONE | www.plosone.org 



6 



November 2013 | Volume 8 | Issue 11 | e81439 



Unexpected Bacterial Pathogens in Ticks 



3. Andcrsson M, Rabcrg L (2011) Wild rodents and novel human pathogen 
candidatus Neoehrliehia mikurensis, Southern Sweden. Emerg Infect Dis 17: 
1716-1718. 

4. Haselbarth K, Tenter AM, Brade V, Krieger G, Hunfcld KP (2007) First ease of 
human babesiosis in Germany - Glinieal presentation and moleeular 
characterisation of the pathogen. Int J Med Microbiol 297: 197-204. 

5. Parola P, Raoult D (2001) Tick-borne bacterial diseases emerging in Europe. 
Glin Microbiol Infect 7: 80-83. 

6. Hubalek Z (2009) Epidemiology of lyme borrcliosis. Curr Probl Dermatol 37: 
31-50. 

7. Piesman J, Eiscn L (2008) Prevention of tick-borne diseases. Annu Rev Entomol 

53: 323-343. 

8. Kik M, Nijhof AM, Balk JA, Jongejan F (2011) Babesia sp. EU1 infection in a 
forest reindeer, The Netherlands. Emerg Infect Dis 17: 936—938. 

9. Mayne PJ (2011) Emerging incidence of Lyme borreliosis, babesiosis, 
bartonellosis, and granulocytic ehrlichiosis in Australia. Int J Gen Med 4: 
845-852. 

10. Fcder HM, Jr., Johnson BJ, O'Connell S, Shapiro ED, Steere AG, et al. (2007) A 
critical appraisal of "chronic Lyme disease". N Engl J Med 357: 1422—1430. 

11. Sauvage V, Gheval J, Foulongne V, Gouilh MA, Paricnte K, et al. (2011) 
Identification of the first human Gyrovirus, a virus related to chicken anemia 
virus. J Virol. 

12. Sauvage V, Foulongne V, Gheval J, Ar Gouilh M, Pariente K, et al. (2011) 
Human polyomavirus related to African green monkey lymphotropic poly- 
omavirus. Emerg Infect Dis 17: 1364—1370. 

13. Cox-Foster DL, Conlan S, Holmes EC, Palacios G, Evans JD, ct al. (2007) A 
metagenomic survey of microbes in honey bee colony collapse disorder. Science 
318:283-287. 

14. Buckwalter MR, Nga PT, Gouilh MA, Fiette L, Bureau JF, et al. (2011) 
Identification of a Novel Ncuropathogcnic Thciler's Murine Encephalomyelitis 
Virus. J Virol 85: 6893-6905. 

1 5. Nakamura S, Maeda N, Miron IM, Yoh M, Izutsu K, et al. (2008) Metagenomic 
diagnosis of bacterial infections. Emerg Infect Dis 14: 1784—1786. 

16. Nakamura S, Nakaya T, Iida T (2011) Metagenomic analysis of bacterial 
infections by means of high-throughput DNA sequencing. Exp Biol Med 
(Maywood) 236: 968-971. 

17. Gheval J, Sauvage V, Frangcul L, Dacheux L, Guigon G, et al. (2011) 
Evaluation of high-throughput sequencing for identifying known and unknown 
viruses in biological samples. J Clin Microbiol 49: 3268-3275. 

18. Hoffmann B, Scheuch M, Hopcr D, Jungblut R, Holsteg M, ct al. (2012) Novel 
orthobunyavirus in cattle, curope, 2011. Emerg Infect Dis 18: 469^472. 

19. Carpi G, Cagnacci F, Wittekindt NE, Zhao F, Qi J, et al. (201 1) Metagenomic 
profile of the bacterial communities associated with Ixodes rieinus ticks. PLoS 
One 6: c25604. 

20. Nakao R, Abe T, Nijhof AM, Yamamoto S, Jongejan F, ct al. (2013) A novel 
approach, based on BLSOMs (Batch Learning Self-Organizing Maps), to the 
microbiome analysis of ticks. ISMEJ 7: 1003-1015. 

21. Reis C, Cote M, Paul RE, Bonnet S (201 1) Questing ticks in suburban forest are 
infected by at least six tick-borne pathogens. Vector Borne Zoonotic Dis 1 1 : 
907-916. 

22. Reis C, Cote M, Le Rhun D, Lecuelle B, Levin ML, et al. (2011) Vector 
competence of the tick Ixodes rieinus for transmission of Bartonella birtlesii. 
PLoS Negl Trop Dis 5: el 186. 

23. Hauck WW (1991) Confidence intervals for scroprevalenec determined from 
pooled sera. Ann Epidemiol 1: 277-281. 

24. Cowling DW, Gardner IA, Johnson WO (1999) Comparison of methods for 
estimation of individual-level prevalence based on pooled samples. Prev Vet 
Med 39: 211-225. 

25. Socolovschi C, Mcdiannikov O, Raoult D, Parola P (2009) Update on tick-borne 
bacterial diseases in Europe. Parasite 16: 259—273. 

26. McKiel JA, Bell EJ, Lackman DB (1967) Rickettsia canada: a new member of 
the typhus group of rickettsiae isolated from Haemaphysalis leporispalustris ticks 
in Canada. Can J Microbiol 13: 503-510. 

27. Brouqui P, Parola P, Fournier PE, Raoult D (2007) Spotted fever rickettsioses in 
southern and eastern Europe. FEMS Immunol Med Microbiol 49: 2-12. 



28. Welindcr-Olsson C, Kjellin E, Vaht K, Jacobsson S, Wcnneras C (2010) First 
case of human "Candidatus Neoehrliehia mikurensis" infection in a febrile 
patient with chronic lymphocytic leukemia. J Clin Microbiol 48: 1956-1959. 

29. Vayssier-Taussat M, Le Rhun D, Buffet JP, Maaoui N, Galan M, et al. (2012) 
Candidatus Neoehrliehia mikurensis in Bank Voles, France. Emerg Infect Dis 
18: 2063-2065. 

30. Doyle CK, Labruna MB, Breitschwerdt EB, Tang YW, Gorstvet RE, et al. 
(2005) Detection of medically important Ehrlichia by quantitative multicolor 
TaqMan real-time polymerase chain reaction of the dsb gene. J Mol Diagn 7: 
504-510. 

31. Chomel BB, Boulouis HJ, Breitschwerdt EB, Kasten RW, Vayssier-Taussat M, 
ct al. (2009) Ecological fitness and strategies of adaptation of Bartonella species 
to their hosts and vectors. Vet Res 40: 29. 

32. Janecek E, Mietze A, Goethe R, Schnieder T, Strube G (2012) Bartonella spp. 
infection rate and B. grahamii in ticks. Emerg Infect Dis 18: 1689-1690. 

33. Dietrich F, Schmidgen T, Maggi RG, Richter D, Matuschka FR, et al. (2010) 
Prevalence of Bartonella henselae and Borrelia burgdorferi sensu lato DNA in 
ixodes rieinus ticks in Europe. Appl Environ Microbiol 76: 1395—1398. 

34. Marsot M, Sigaud M, Chapuis JL, Fcrquel E, Cornet M, et al. (201 1) Introduced 
Siberian chipmunks (Tamias sibiricus barberi) harbor more-diverse Borrelia 
burgdorferi sensu lato genospecies than native bank voles (Myodes glareolus). 
Appl Environ Microbiol 77: 5716-5721. 

35. Platonov AE, Karan LS, Kolyasnikova NM, Makhneva NA, Toporkova MG, et 
al. (2011) Humans infected with relapsing fever spirochete Borrelia miyamotoi, 
Russia. Emerg Infect Dis 17: 1816-1823. 

36. Geller J, Nazarova L, Katargina O, Jarvckulg L, Fomenko N, et al. (2012) 
Detection and genetic characterization of relapsing fever spirochete Borrelia 
miyamotoi in Estonian ticks. PLoS One 7: c51914. 

37. Logares R, Sunagawa S, Salazar G, Cornejo-Castillo FM, Ferrera I, ct al. (2013) 
Metagenomic 16S rDNA Illumina tag's are a powerful alternative to amplieon 
sequencing to explore diversity and structure of microbial communities. Environ 
Microbiol. 

38. Li H, Jiang J, Tang F, Sun Y, Li Z, et al. (2013) Wide distribution and genetic 
diversity of "Candidatus Neoehrliehia mikurensis" in rodents from China. Appl 
Environ Microbiol 79: 1024-1027. 

39. Li H, Jiang JF, Liu W, Zheng YC, Huo QB, ct al. (2012) Human infection with 
Candidatus Neoehrliehia mikurensis, China. Emerg Infect Dis 18: 1636—1639. 

40. Fukunaga M, Korcki Y (1995) The flagellin gene of Borrelia miyamotoi sp. nov. 
and its phylogenetic relationship among Borrelia species. FEMS Microbiol Lett 
134: 255-258. 

41. Fukunaga M, Takahashi Y, Tsuruta Y, Matsushita O, Ralph D, et al. (1995) 
Genetic and phenotypic analysis of Borrelia miyamotoi sp. nov., isolated from 
the ixodid tick Ixodes pcrsulcatus, the vector for Lyme disease in Japan. Int J Syst 
Baetcriol 45: 804-810. 

42. Fraenkcl GJ, Garpmo U, Berglund J (2002) Determination of novel Borrelia 
genospecies in Swedish Ixodes rieinus ticks. J Clin Microbiol 40: 3308-3312. 

43. Richter D, Schlee DB, Matuschka FR (2003) Relapsing fever-like spirochetes 
infecting European vector tick of Lyme disease agent. Emerg Infect Dis 9: 697- 
701. 

44. Scoles GA, Papero M, Beati L, Fish D (2001) A relapsing fever group spirochete 
transmitted by Ixodes scapularis ticks. Vector Borne Zoonotic Dis 1: 21-34. 

45. Sarksyan DS, Platonov AE, Karan LS, Malinin IE, Khalitova LI, et al. (2012) 
[Clinical presentation of "new" tick-borne borrcliosis caused by Borrelia 
miyamotoi] . Ter Arkh 84: 34-^-1. 

46. Cotte V, Bonnet S, Le Rhun D, Le Naour E, Ghauvin A, ct al. (2008) 
Transmission of Bartonella henselae by Ixodes rieinus. Emerg Infect Dis 14: 
1074-1080. 

47. Klee SR, TyczkaJ, Ellerbrok H, Franz T, Linke S, et al. (2006) Highly sensitive 
real-time PGR for specific detection and quantification of Coxiella burnetii. 
BMC Microbiol 6: 2. 

48. Vcrsage JL, Severin DD, Chu MC, Petersen JM (2003) Development of a 
multitarget real-time TaqMan PGR assay for enhanced detection of Franeisella 
tularensis in complex specimens. J Clin Microbiol 41: 5492—5499. 

49. Socolovschi C, Pages F, Raoult D (2012) Rickettsia felis in Aedes albopictus 
mosquitoes, Libreville, Gabon. Emerg Infect Dis 18: 1687-1689. 



PLOS ONE | www.plosone.org 



7 



November 2013 | Volume 8 | Issue 11 | e81439 



